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Abstract 

The search for experimental demonstrations of the quantum behavior of macroscopic 
mechanical resonators is a fastly growing field of investigation and recent results 
suggest that the generation of quantum states of resonators with a mass at the 
microgram scale is within reach. In this chapter we give an overview of two important 
topics within this research field: cooling to the motional ground state, and the 
generation of entanglement involving mechanical, optical and atomic degrees of 
freedom. We focus on optomechanical systems where the resonator is coupled to 
one or more driven cavity modes by the radiation pressure interaction. We show 
that robust stationary entanglement between the mechanical resonator and the 
output fields of the cavity can be generated, and that this entanglement can be 
transferred to atomic ensembles placed within the cavity. These results show that 
optomechanical devices are interesting candidates for the realization of quantum 
memories and interfaces for continuous variable quantum communication networks. 
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1 Introduction 

Mechanical resonators at the micro- and nano-meter scale are widely employed 
for a large variety of applications, more commonly as sensors or actuators in 
integrated electrical, optical, and opto-electronical systems p^|l2|l3|l^ . Modifica- 
tions of the resonator motion can be detected with high sensitivity by looking 
at the radiation (or electric current) which interacted with the resonator. For 
example, small masses can be detected by measuring the frequency shift in- 
duced on the resonator, while tiny displacements (or weak forces inducing such 
displacements) can be measured by detecting the corresponding phase shift of 
the light interacting with it [2]. The resonators are always subject to thermal 
noise, which is due to the coupling with internal and/or external degrees of 
freedom and is one of the main factors limiting the sensitivity of these devices. 
However, due to the progress in nanofabrication techniques, the mechanical 
quality factor Qm (which quantifies this undesired coupling to environmental 
degrees of freedom) is steadily improving, suggesting that in the near future 
these devices will reach the regime in which their sensitivity is limited by the 
ultimate quantum limits set by the Heisenberg principle. The importance of 
the limits imposed by quantum mechanics on the resonator motion was first 
pointed out by Braginsky and coworkers [5] in the completely different con- 
text of massive resonators employed in the detection of gravitational waves 
[6]. However, in recent years the quest for the experimental demonstration 
of genuine quantum states of macroscopic mechanical resonators has spread 
well beyond the gravitational wave physics community and has attracted a 
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wide interest. In fact, the detection of an unambiguous signature of the quan- 
tum behavior of a macroscopic oscillator, with a mass at least of the order of 
a microgram, would shed further light onto the quantum-classical boundary 
problem [7]. In fact, nothing in the principles of quantum mechanics prevents 
macroscopic systems to be prepared in genuine quantum states. However, it 
is not yet clear how far one can go in this direction [8], and a complete un- 
derstanding of how classical behavior emerges from the quantum substrate re- 
quires the design and the implementation of dedicated experiments. Examples 
of this kind are single-particle interference of macro- molecules [5], the demon- 
stration of entanglement between collective spins of atomic ensembles [TU] . 
and of entanglement in Josephson-j unction qubits [llj. For what concerns me- 
chanical resonators, the experimental efforts are currently focusing on cooling 
them down to their motional ground state [2J. This goal has not been achieved 
yet, but promising results in this direction have been obtained in different 

different examples of mechanical resonators coupled either to radiative or to 
electrical degrees of freedom. Ground state cooling of microgram-scale res- 
onators seems to be within reach, as already suggested by various theoretical 

ical oscillator can be coupled to another system so that the latter can act 
as an effective zero-temperature reservoir. In the first part of this chapter 
we shall review the problem of ground state cooling of a mechanical res- 
onator, by focusing onto the case where the role of effective zero-temperature 
"fridge" is played by an optical cavity mode, coupled to the resonator by ra- 
diation pressure. In this case this interaction can be exploited for cooling in 
two different ways: i) back-action, or self-cooling [33ll39|40|41|42|43j in which 
the off-resonant operation of the cavity results in a retarded back action on 
the mechanical system and hence in a "self" -modification of its dynamics 

where the oscillator position is measured through a phase-sensitive detection 
of the cavity output and the resulting photocurrent is used for a real-time 
correction of the dynamics P^16lll9f22ll28j . We shall compare the two ap- 
proaches and see that while back-action cooling is optimized in the good cavity 
limit where the resonator frequency is larger than the cavity bandwidth, cold 
damping is preferable in the opposite regime of larger cavity bandwidths [H] . 
It should be noticed that the model Hamiltonian based on radiation pres- 
sure coupling between an optical cavity mode and one movable cavity mirror 
is quite general and immediately extendible to other situations, such as the 
toroidal microcavities of Refs. [20|25] . the capacitively coupled systems of 
Refs. and even atomic condensate systems [H]. 

From the theory side, the generation of other examples of quantum states of 
a micro-mechanical resonator has been also considered. The most relevant ex- 
amples are given by squeezed and resonator-field (or atoms) entangled states. 
Squeezed states of nano-mechanical resonators are potentially useful for 
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surpassing the standard quantum limit for position and force detection [5], and 
could be generated in different ways, either by coupling to a qubit [16], or by 
measurement and feedback schemes [SBfTT] . Entanglement is instead the char- 
acteristic element of quantum theory, because it is responsible for correlations 
between observables that cannot be understood on the basis of local realistic 
theories [48J. For this reason, there has been an increasing interest in establish- 
ing the conditions under which entanglement between macroscopic objects can 
arise. Relevant experimental demonstration in this directions are given by the 
entanglement between collective spins of atomic ensembles (TUI, and between 
Josephson-junction qubits [TT]. Then, starting from the proposal of Ref. [IH] in 
which two mirrors of a ring cavity are entangled by the radiation pressure of the 
cavity mode, many proposals involved nano- and micro-mechanical resonators, 
eventually entangled with other systems. One could entangle a nanomechanical 
oscillator with a Cooper-pair box [50j, while Ref. ^5]J studied how to entangle 
an array of nanomechanical oscillators. Further proposals suggested to entan- 
gle two charge qubits [52] or two Josephson junctions [5^ via nanomechanical 
resonators, or to entangle two nanomechanical resonators via trapped ions 
[51], Cooper pair boxes [55], or dc-SQUIDS [56]. More recently, schemes for 
entangling a superconducting coplanar waveguide field with a nanomechanical 
resonator, either via a Cooper pair box within the waveguide [27], or via direct 
capacitive coupling [58], have been proposed. After Ref. 09], other optome- 
chanical systems have been proposed for entangling optical and/or mechanical 
modes by means of the radiation pressure interaction. Ref. [59] considered two 
mirrors of two different cavities illuminated with entangled light beams, while 
Refs. [501I61II62|[S5] considered different examples of double-cavity systems in 
which entanglement either between different mechanical modes, or between 
a cavity mode and a vibrational mode of a cavity mirror have been studied. 
Refs. [M1I65] considered the simplest scheme capable of generating stationary 
optomechanical entanglement, i.e., a single Fabry- Perot cavity either with one 
[64], or both [HS], movable mirrors. 

In the second part of the chapter we shall focus on the generation of stationary 
entanglement by starting from the Fabry- Perot model of Ref. [64], which is re- 
markable for its simplicity and robustness against temperature, and extend its 
study in various directions. In fact, entangled optomechanical systems could 
be profitably used for the realization of quantum communication networks, 
in which the mechanical modes play the role of local nodes where quantum 
information can be stored and retrieved, and optical modes carry this informa- 
tion between the nodes. Refs. [66ll67|68j proposed a scheme of this kind, based 
on free-space light modes scattered by a single reflecting mirror, which could 
allow the implementation of continuous variable (CV) quantum teleportation 
[66], quantum telecloning [67], and entanglement swapping [68]. Therefore, any 
quantum communication application involves traveling output modes rather 
than intracavity ones, and it is important to study how the optomechanical 
entanglement generated within the cavity is transferred to the output field. 
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Furthermore, by considering the output field, one can adopt a multiplexing 
approach because, by means of spectral filters, one can always select many dif- 
ferent traveling output modes originating from a single intracavity mode. One 
can therefore manipulate a multipartite system, eventually possessing mul- 
tipartite entanglement. We shall develop a general theory showing how the 
entanglement between the mechanical resonator and optical output modes 
can be properly defined and calculated [69]. We shall see that, together with 
its output field, the single Fabry-Perot cavity system of Ref. |64j represents 
the "cavity version" of the free-space scheme of Refs. [BSEZ!- In fact, as it 
happens in this latter scheme, all the relevant dynamics induced by radia- 
tion pressure interaction is carried by the two output modes corresponding 
to the first Stokes and anti-Stokes sidebands of the driving laser. In particu- 
lar, the optomechanical entanglement with the intracavity mode is optimally 
transferred to the output Stokes sideband mode, which is however robustly en- 
tangled also with the anti-Stokes output mode. We shall see that the present 
Fabry-Perot cavity system is preferable with respect to the free space model 
of Refs. [66|67] . because entanglement is achievable in a much more accessi- 
ble experimental parameter region. We shall then extend the analysis to the 
case of a doubly-driven cavity mode. We shall see that a peculiar parameter 
regime exists where the optomechanical system, owing to the combined ac- 
tion of the two driven modes, is always stable and is characterized by robust 
entanglement between the resonator and the cavity output fields. 

In the last Section we shall investigate the possibility to couple and entangle in 
a robust way optomechanical systems to atomic ensembles, in order to achieve 
a strongly-coupled hybrid multipartite system [7DH7T] . We shall see that this 
is indeed possible, especially when the atomic ensemble is resonant with the 
Stokes sideband induced by the resonator motion. Such hybrid systems might 
represent an important candidate for the realization of CV quantum interfaces 
within CV quantum information networks. 



2 Cavity optomechanics via radiation pressure 

The simplest cavity optomechanical system consists of a Fabry-Perot cavity 
with one heavy, fixed mirror through which a laser of frequency ui drives a 
cavity mode, and another light end-mirror of mass m (typically in the micro 
or nanogram range), free to oscillate at some mechanical frequency Um- Our 
treatment is however valid also for other cavity geometries in which one has an 
optical mode coupled by radiation pressure to a mechanical degree of freedom. 
A notable example is provided by silica toroidal optical microcavities which 
are coupled to radial vibrational modes of the supporting structure [20|72] . 
Radiation pressure typically excites several mechanical degrees of freedom of 
the system with different resonant frequencies. However, a single mechanical 
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mode can be considered when a bandpass filter in the detection scheme is used 
[73] and couphng between the different vibrational modes can be neglected. 
One has to consider more than one mechanical mode only when two close 
mechanical resonances fall within the detection bandwidth (see Ref. j74j for 
the effect of a nearby mechanical mode on cooling and entanglement). The 
Hamiltonian of the system describes two harmonic oscillators coupled via the 
radiation pressure interaction, and reads [75] 



H = huj^a)a + + q^) - HGoa^aq + ikSia^e-''^'' - ae'"^''). (1) 



The first term describes the energy of the cavity mode, with lowering operator 
a {[a, = 1), frequency Uc (and therefore detuned by Aq = Uc — oji from the 
laser), and decay rate k. The second term gives the energy of the mechanical 
mode, described by dimensionless position and momentum operators q and p 
{[liP] = '^)- The third term is the radiation-pressure coupling of rate Gq = 



{ujc/ L)^Jh/muJm, where m is the effective mass of the mechanical mode [73j, 
and L is an effective length that depends upon the cavity geometry: it coincides 
with the cavity length in the Fabry-Perot case, and with the toroid radius 
in the case of Refs. [20|I72] . The last term describes the input driving by a 
laser with frequency ui, where S is related to the input laser power V by 



\S\ = y2VK/huJi. One can adopt the single cavity mode description of Eq. ([T]) 
as long as one drives only one cavity mode and the mechanical frequency Um 
is much smaller than the cavity free spectral range FSR ~ c/2L. In this case, 
in fact, scattering of photons from the driven mode into other cavity modes 
is negligible [76] . 



2. 1 Langevin equations formalism 

The djTiamics are also determined by the fluctuation-dissipation processes af- 
fecting both the optical and the mechanical mode. They can be taken into 
account in a fully consistent way [75] by considering the following set of non- 
linear QLE (quantum Langevin equations), written in a frame rotating at 

q = uJmP, (2) 
p = -ujmq - imP + Goo^a + ^, (3) 
a = -{K + iAQ)a + iGQaq + £ + V^j^a^"". (4) 

The mechanical mode is affected by a viscous force with damping rate 7^, 
and by a Brownian stochastic force with zero mean value i{t), possessing the 



1 



2 




6 



correlation function [75|77] 



mm 



7^ f ^^-^^it-t')^ 



2ti 



coth 



fUjJ 



2kBTo, 



(5) 



where ks is the Bohzmann constant and Tq is the temperature of the reservoir 
of the micromechanical oscillator. The correlation function and the commuta- 
tor of the Gaussian stochastic force ^ (t) are not proportional to a Dirac delta 
and therefore ^{t) is a non-Markovian stochastic process. This fact guaran- 
tees that the QLE of Eqs. ©-(Ill) preserve the correct commutation relations 
between operators during the time evolution [75]. However, a Markovian de- 
scription of the symmetrized correlations of ^(t) is justified in two different 
limits, which are both met in typical experimental situations: i) not too low 
temperatures kBTo/huJm ^ 1, which for typical values is satisfied even at cryo- 
genic temperatures; ii) high mechanical quality factor Q = Um/'Jm oo [78] . 
which is an important condition for the observation of quantum effects on the 
mechanical resonator. In this case the correlation function of Eq. ([5D can be 
approximated as 



mm) 



In 



(2no + l)6{t -t')+t 



S'it-t') 



(6) 



where uq = {exp{fiuJm/kBTo} — 1) is the mean thermal excitation number 
of the resonator and 6'{t — t') denotes the derivative of the Dirac delta. 



The cavity mode amplitude instead decays at the rate k and is affected by the 
vacuum radiation input noise a*"(t), whose correlation functions are given by 
179] 



{d-{t)a'-'\t')) = [Nico,) + 1] 5it - t'). (7) 
{a'^'^{t)a'^t'))=N{uj,)6{t-t'), (8) 

where N{ujc) = (expl^c/^s^o} ~ is the equilibrium mean thermal pho- 
ton number. At optical frequencies hwc/kBTo ^ 1 and therefore N{ujc) — 0, 
so that only the correlation function of Eq. ([7j) is relevant. 

Equations ©-(jl]) are not easy to analyze owing to the nonlinearity. However, 
one can proceed with a linearization of operators around the steady state. The 
semiclassical steady state is characterized by an intracavity field amplitude as 
{\as\ S> 1), and a new equilibrium position for the oscillator, displaced by g^- 
The parameters as and Qs are the solutions of the nonlinear algebraic equations 
obtained by factorizing Eqs. ©-(jlj) and setting the time derivatives to zero: 
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(9) 



where the latter equation is in fact the nonhnear equation determining a^, 
since the effective cavity detuning A, including radiation pressure effects, is 
given by [80] 

r^n, |2 

A = Ao-^^^^. (11) 

Rewriting each Heisenberg operator of Eqs. ([2])-(jl]) as the c-number steady 
state value plus an additional fluctuation operator with zero mean value, one 
gets the exact QLE for the fluctuations 



Sq = uJmSp, (12) 
Sp = —uJmSq — '~)m^P + Co {cisSa) + a*(5a^ + 6a) 8a + ^, (13) 
5a = -{k + iA)6a + iGo (a^ + Sa) 6q + v^a*". (14) 



Since we have assumed la^l ^ 1, one can safely neglect the nonlinear terms 
5a^(5a and SaSq in the equations above, and get the linearized QLE 



5q = uj^5p, (15) 

6p = -uJmSq - 'jmSp + G6X + ^, (16) 

6X = -k5X + A6Y + v^X™, (17) 

SY = -k5Y - A5X + G5q + ^PIhY' 



Here we have chosen the phase reference of the cavity field so that as is real and 
positive, we have defined the cavity field quadratures bX = (Sa + 5a^)/^/2 and 
5Y = (Sa — 5a^)/i\/2 and the corresponding Hermitian input noise operators 
Xin ^ (^m ^ am,t^/^ yin ^ ^^in _ a^n^y-^^ rj^^ linearized QLE 

show that the mechanical mode is coupled to the cavity mode quadrature 
fluctuations by the effective optomechanical coupling 



which can be made very large by increasing the intracavity amplitude a^- 
Notice that together with the condition uJm c/L which is required for the 
single cavity mode description, la^l ^ 1 is the only assumption required by 
the linearized approach. This is in contrast with the perturbative approaches 
described in [30] , where a reduced master equation of the mechanical resonator 
is derived under the weak-coupling assumption G <^ ujrn- 
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2.2 Stability analysis 



The stability analysis can be performed on the linearized set of equations 
Eqs. ©-(jlD by using the Routh-Hurwitz criterion [ST]. Two conditions are 
obtained 



Sl = 27mK I + {uJm - A)^ + {uJm + A) 



+ 7r7 



(20) 

(7„ + 2k) [k^ + A^) + 2kujI] } + Auj.^G^ (7^ + 2k)' > 0, (21) 

(22) 



S2 = (/t' + A') - G'^A > 0. 



The violation of the first condition, si < 0, indicates instability in the domain 
of blue-detuned laser (A < 0) and it corresponds to the emergence of a self- 
sustained oscillation regime where the mirror effective damping rate vanishes. 
In this regime, the laser field energy leaks into field harmonics at frequencies 
uJi±ruJm {r = 1, 2...) and also feeds the mirror coherent oscillations. A complex 
multistable regime can emerge as described in [82] . The violation of the second 
condition S2 < indicates the emergence of the well-known effect of bistable 
behavior observed in [H3] and occurs only for positive detunings (A > 0). In 
the following we restrict our analysis to positive detunings in the stable regime 
where both si and S2 conditions are fulfilled. A parametric plot showing the 
domain of stability in the red-detuning regime A > is shown in Fig. [1] where 
we have plotted the stability parameter 



V 



UJm («:' + A2) 



(23) 



Negative values of r] indicate the emergence of instability. We have chosen 
the following set of parameters which will be used extensively throughout the 



chapter and which is denoted by po= 



.Qm,m,L,X,,To) = (27r x 10 MHz, 



10^, 30 ng, 0.5 mm, 1064 nm, 0.6 K). These values are comparable to those 
used in recent experiments 



2.3 Covariance matrix and logarithmic negativity 



The mechanical and intracavity optical mode form a bipartite continuous vari- 
able (CV) system. We are interested in the properties of its steady state which, 
due to the linearized treatment and to the Gaussian nature of the noise op- 
erators, is a zero-mean Gaussian state, completely characterized by its sym- 
metrized covariance matrix (CM). The latter is given by the 4x4 matrix with 
elements 

{Ui {00) Um {00) +Um {00) Ui {00)) 

yim = , l^4j 
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Fig. 1. Stability condition in the red-detuning region, (a) Contour plot of the stabil- 
ity parameter r] of Eq. (j23p as a function of input power V and normalized detuning 
/S./uJm- The parameter set pQ={u;m,Qm,i^, L, Xc,Tq) = {2it x 10 MHz, 10^, 30 ng, 
0.5 mm, 1064 nm, 0.6 K has been used, together with = 8 x 10^ (corresponding 
to K = O.STwm). The blue area corresponds to the unstable regime, (b) Stability 
parameter rj versus V and the normalized cavity decay rate K,/ujm at A = 

where Um{oo) is the asymptotic value of the m-th component of the vector of 
quadrature fluctuations 



u{t) = {5q{t),5p{t),5X{t),5Y{t)y . 



(25) 



Its time evolution is given by Eqs. ( ]T5l) - (]T8ll . which can be rewritten in compact 
form as 

'^-u{t)= Au{t) + v{t), (26) 



dt 



with A the drift matrix 



A 
















'~1m 


G 











— K 


A 


G 





-A 


— K 



(27) 



and v{t) the vector of noises 

v{t) = (0,e(t),v^X^"(t),v^r"(t))\ (28) 
The steady state CM can be determined by solving the Lyapunov equation 

AV + VAT = (29) 
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where D is the 4x4 diffusion matrix which characterizes the noise correlations 
and is defined by the relation {rii (t) rim (t') + rim (t') ni(t)) /2 = Dimwit — t'). 
Using Eqs. ([H])-(IZD, D can be written as 

D = diag[0,7m(2no + !),«:,«:]. (30) 

Eq. (12^ is a linear equation for V and it can be straightforwardly solved, but 
the general exact expression is very cumbersome and will not be reported here. 

The CM allows to calculate also the entanglement of the steady state. We 
adopt as entanglement measure the logarithmic negativity Ej^, which is de- 
fined as [84j 

Eat = max[0,-ln2r/"]. (31) 

1 /2 

Here r]- = 2'^/^ [s(V) - [S(V)2 - 4 det V]^/^] and S(V) = detVi + detVs- 
2det Vc, with Vi, V2 and Vc being 2x2 block matrices of 



V 



(32) 



A bimodal Gaussian state is entangled if and only if rj^ < 1/2, which is 
equivalent to Simon's necessary and sufficient entanglement non-positive par- 
tial transpose criterion for Gaussian states [85], which can be written as 
4detV < S(V) — 1/4. Logarithmic negativity is a convenient entanglement 
measure because it is the only one which can always be explicitly computed 
and it is also additive. The drawback of Ej^ is that, differently from the en- 
tanglement of formation and the distillable entanglement, it is not strongly 
super-additive and therefore it cannot be used to provide lower-bound esti- 
mates of the entanglement of a generic state by evaluating the entanglement 
of Gaussian state with the same correlation matrix This fact however is 
not important in our case because the steady state of the system is Gaussian 
within the validity limit of our linearization procedure. 



3 Ground state cooling 



The steady state CM V determines also the mean energy of the mechanical 
resonator, which is given by 



U 



fku„ 



[Vii + V22] = n + 



(33) 



where n = {exp{hujm/ ksT} — 1)~ is the occupancy corresponding to a bath 
temperature T. Obviously, in the absence of coupling to the cavity field it is 
n = riQ, where riQ corresponds to the actual temperature of the environment Tq. 
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Fig. 2. Setup for feedback cooling (cold damping). The cavity output field is homo- 
dyne detected (thus acquiring information about the mirror position) and a force 
proportional to its derivative is fed back to the mirror. 

The optomechanical coupling with the cavity mode can be used to 'engineer' 
an effective bath of much lower temperature T -C Tq, so that the mechanical 
resonator is cooled. Let us see when it is possible to reach the ideal condition 
n -C 1, which corresponds to ground state cooling. 



3.1 Feedback cooling 

A simple way for cooling an object is to continuously detect its momentum 
and apply 'corrective kicks' that continuously reduce it eventually to zero 
[32|35|36j . This is the idea of feedback cooling illustrated in Fig. [2] where 
the mirror position is detected via phase-sensitive homodyne detection of the 
cavity output field and a force proportional to the time derivative of the output 
signal (thus to the velocity) is fed back to it. By Fourier transforming Eq. ( |T8i) 
one obtains 



which shows that the intracavity phase-quadrature is sensitive to the mirror 
motion and moreover its optimal sensitivity is reached at resonance, when 
A = 0. In this latter condition 5X{uj) is not sensitive to the mirror motion, 
suggesting that the strongest feedback effect is obtained by detecting the out- 
put phase-quadrature and feeding it back to the resonator. 



G(k — iuj) 



Sq{uj) + noise terms. 



(34) 
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3.1.1 Phase quadrature feedback 

As a consequence we set A = and add a feedback force in Eq. (fT6|) so that 

5p = -uJmSq - imh + G5X + i- f dsg{t - s)6Y^'\s), (35) 

J —oo 

where Y^^^{s) is the estimated intracavity phase-quadrature, which, using 
input-output relations [79] and focusing onto the ideal scenario of perfect 
detection, is given by 

SY''\t) = = 6Y{t) - (36) 



/2k V2K 

The filter function g{t) is a causal kernel and g{uj) is its Fourier transform. 
We choose a simple standard derivative high-pass filter 



d 



lujgcd 



9{t) = gcaj^ e{t)ojf,e'-f^' g{to) = (37) 



1 — iio/ujfb^ 



so that ujji^ plays the role of the time delay of the feedback loop, and gcd > 
is the feedback gain. The ideal derivative limit is obtained when Ufb oo, 
implying g^u) = —iojgcd and therefore g{t) = gcd^'if). In this limit the feedback 
force is equal (apart from an additional noise term) to —gcdSY which, due to 
Eq. fl51|) . is an additional viscous force —{gcdG/K)6q only in the bad cavity 
limit K > 

One can solve the Langevin equations supplemented with the feedback term 
in the Fourier domain. In fact, the two steady state oscillator variances {Sq"^) 
and {Sp"^) can be expressed by the following frequency integrals 

(38, 

where Sg'^{u) is the position noise spectrum. Its explicit expression is given by 

Sf{^) = |Xe?/(^)r[5*/.(^) + SrAoo) + Sfbioo)], (39) 

where the thermal, radiation pressure and feedback-induced contributions are 
respectively given by 



Sth{uj) = ^^^^ coth 



m 



' tVjJ 



Q'. 

Srp{w) = 



K 



2 



Q / N 1^(^)1 



(40) 
(41) 
(42) 
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and Xeffi^) is the susceptibility of the mechanical oscillator modified by the 
feedback 



-Mu) ul-u'- .a;7™ + ^M^j ' . (43) 

K — tUJ 



This effective susceptibility contains the relevant physics of cold damping. In 
fact it can be rewritten as the susceptibility of an harmonic oscillator with 
effective (frequency-dependent) damping and oscillation frequency. The modi- 
fication of resonance frequency (optical spring effect [21ii36j ) is typically small 
for the chosen parameter regime — 1 MHz) and the only relevant effect 
of feedback is the modification of the mechanical damping which, in the case 
of the choice of Eq. (I57|) , is given by 

This expression shows that the damping of the oscillator may be significantly 
increased due to the combined action of feedback and of radiation pressure 
coupling to the field. In the ideal limit of instantaneous feedback and of a bad 
cavity, k, tu/b ^ Um, 7™,, effective damping is frequency-independent and given 
by j^'^''^'^ — 7m + QcdG^ml = 7m(l + 92), where we have defined the scaled, 
dimensionless feedback gain g2 = gcdGuJml i^lm [36] . 

The presence of cold-damping feedback also modifies the stability conditions. 
The Routh-Hurwitz criteria are equivalent to the conditions that all the poles 
of the effective susceptibility of Eq. ( I43l) are in the lower complex half-plane. 
For the choice of Eq. (l37j) there is only one non-trivial stability condition, 
which reads 



Scd = 



^mKUJfb + gcdGuJmUJfb + + ^ fb) [(/« + 7m) (/^ + t^/fe)(7m + ^ fb) 

\ - g^dGuJrn^^fb] - Kul^UJfbiK + 7^ + UJ fb^ > 0. (45) 



This condition shows that the system may become unstable for large gain 
and finite feedback delay-time and cavity bandwidth because in this limit the 
feedback force can be out-of-phase with the oscillator motion and become an 
accelerating rather than a viscous force 



The performance of cold-damping feedback for reaching ground state cooling 
is analyzed in detail in Ref. [41], which shows that the optimal parameter 
regime is k t^/b ~ Um ^ 7m, which correspond to a bad-cavity limit and a 
finite-bandwidth feedback, i.e., with a feedback delay-time comparable to the 
resonator frequency. One gets in this case 
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(46) 



(47) 



where we have defined the scaled dimensionless input power ( = 2G'^/K'jrn- 
These two expressions show that with cold-damping feedback, ^ 
i.e., energy equipartition does not hold anymore. The best cooling regime is 
achieved for Ufb ~ SuJm and — ^ (i.e. Qcd — '^G/uim), i-e. for large hut 
finite feedback gain [35]l36f41] . This is consistent with the fact that stability 
imposes an upper bound to the feedback gain when k and ujfb are finite. 
The optimal cooling regime for cold damping is illustrated in Fig. [3^, where 
n is plotted versus the feedback gain g^d and the input power P, at fixed 
K = Scjm (bad-cavity condition) and cj/f, = S.ScJm- Fig. [3b instead explicitly 
shows the violation of the equipartition condition even in this regime close to 
ground state (the feedback gain is fixed at the value gcd = 1-2): the resonator 
is in a position-squeezed thermal state corresponding to a very low effective 
temperature. 



3.1.2 Generalized quadrature feedback 

The above analysis shows that cold-damping feedback better cools the me- 
chanical resonator when the feedback is not instantaneous and therefore the 
feedback force is not a simple viscous force. This suggests that one can further 
optimize feedback cooling by considering a generalized estimated quadrature 
which is a combination of phase and amplitude field quadratures. In fact one 
may expect that in the optimal regime, the information provided by the am- 
plitude quadrature X°"*(t) is also useful. 

Therefore, in order to optimize cooling via feedback, we apply a feedback force 
involving a generalized estimated quadrature 

^yest^^^ ^ Y"^\t) COS 6 + X"-\t) Sin 9 ^ ^^^^ 

which is a linear combination of y°"*(t) and X"^^{t) and where ^ is a detection 
phase which has to be optimized. The adoption of the new estimated quadra- 
ture leads to three effects: i) a modification of the expression for Xe'ffi^) of 
Eq. (Il3ll where g{uj) is replaced by g^uj) cos9; ii) a consequent reduction of the 
feedback-induced shot noise term Sfb{uj); iii) a reduction of radiation pressure 
noise. In fact, the radiation pressure and feedback- induced noise contributions 
become 
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Fig. 3. Feedback cooling, (a) Contour plot of function of V and gcd- The 

parameters are po, k = 5uJm and ujfb = S.Swm- (b) Illustration of the violation of 
energy equipartition around the optimal cooling regime. Parameters as before with 
Qcd = 1-2. (c) n versus the phase of the generalized quadrature 9 for two sets of gcd 
and V: the (upper) blue curve corresponds to Qcd = 0.8 and "P = 20 mW, while the 
(lower) red curve corresponds to gcd = 1-2 and V = 50 mW. (d) Comparison of n 
versus the input power V between the case of standard cold damping feedback ^ = 
(upper red curve) and at a generalized detected quadrature with phase 6 = O.IStt 
(lower blue curve). Parameters as before, with gcd = 1-2. 

2 



G'k 



1 - 



g{uj) sin 6* 



2GKi 



K + iu] 



\9{^)\' 
4k 



cos^ e. 



(49) 
(50) 



An improvement over the standard cold-damping feedback scheme can be 
obtained when the shot noise reduction effect predominates over the reduction 
of the effective damping due to feedback. This can be seen in Fig. (jSt) where 
for two different choices for gcd and V, the occupancy n is plotted versus 9. For 
one of these optimal phases, 9opt = O.lSvr, we plot in Fig. (jSji) n as a function of 
V and compare it with the results of the standard phase quadrature feedback 
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to conclude that improvement via detection of a rotated output quadrature is 
indeed possible. 



3.2 Back-action cooling 



In analogy with well-known methods of atom and ion cooling [87f88] . one 
can also think of cooling the mechanical resonator by exploiting its coherent 
coupling to a fast decaying system which provides an additional dissipation 
channel and thus cooling. In the present situation, radiation pressure couples 
the resonator with the cavity mode and the fast decaying channel is provided 
by the cavity photon loss rate k. An equivalent description of the process can 
be given in terms of dynamical backaction : the cavity reacts with a delay 
to the mirror motion and induces correlations between the radiation pressure 
force and the Brownian motion that lead to cooling or amplification, depending 
on the laser detuning. A quantitative description is provided by considering 
scattering of laser photons into the motional sidebands induced by the mirror 
motion (see Fig. H]) [S^4Uf41j . Stokes (red) and anti-Stokes (blue) sidebands 
are generated in the cavity at frequencies uji±ujrn- Laser photons are scattered 
by the moving oscillator into the two sidebands with rates 

A+ = — (51) 



+ (A ± UJmf 

simultaneously with the absorption (Stokes, A^) or emission (anti-Stokes, A^) 
of vibrational phonons. The inequality A- > ^4+ leads to a decrease in the 
oscillator phonon occupation number and thus to cooling. Eq. flSTl) shows that 
this occurs when A > and that an effective optical cooling rate, 

T = A -A = 2G'Au;^K 

" + [k2 + (^^-A)2][/c2 + (^^ + A)2]' ^ ' 

can be defined, providing a measure of the coupling rate of the resonator 
with the effective zero-temperature environment represented by the decaying 
cavity mode. Since the mechanical damping rate 7^ is the coupling rate with 
the thermal reservoir of the resonator, one can already estimate that, when 
r ^ 7m, the mechanical oscillator is cooled at the new temperature T 
(7„/r) To. 

One can perform a more precise and rigorous derivation of the cooling rate 
and steady state occupancy by using Eq. fl5^ . The position and momentum 
variances can be in fact obtained by solving Eq. (!29|) or, equivalently, by 
solving the linearized QLE in the Fourier domain and integrating the resulting 
noise spectra. The result of these calculations, in the limit of large mechanical 
quality factor Qm, reads 
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Fig. 4. Setup for cavity backaction cooling. Optical sidebands are scattered unevenly 
by the moving mirror. When the anti-Stokes sideband is resonant with the cavity 
(A = LOm), an effective flow of energy from the mirror out of the cavity leads to an 
effective cooling. 
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2 (A^ - K^) - ul 



A2 



(55) 
(56) 



In the perturbative limit Um ^ nojrn,G and k 3> ■ym,G, Eqs. (I53l)-fl54l) sim- 
plify to {6q^) =i {6p^) ~ n + 1/2, with n ~ [y-mHo + A+] / [7^ + T], which 
reproduces the result of [39P0] . This indicates that ground state cooling is 
reachable when 'jmiT'O < T and provided that the radiation pressure noise con- 
tribution A^/T ~ K^/ (4a;^) is also small. The optical damping rate T can 
be increased by cranking up the input cavity power and thus G. However, 
when one considers the limitations imposed by the stability condition rj > 0, 
one finds that there is an upper bound for G and consequently F. This is 
shown in Figs. [S^-lSt, where one sees that for the chosen parameter regime 
Po, optimal cooling is achieved for A ~ Um (when the anti-Stokes sideband is 
resonant with the cavity, as expected), and in a moderate good-cavity condi- 
tion, n/uJm ~ 0.2. Fig. shows that close to this optimal cooling condition, 
equipartition is soon violated when the input power (and therefore the effective 
coupling G) is further increased: the position variance becomes much larger 
than the momentum variance and it is divergent at the bistability threshold 
(see Eq. (|54l)). 
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Fig. 5. Back action cooling, (a) Contour plot of n versus A/ujm and V. The parame- 
ters are po and k = O.SltOrn- Optimal cooling is seen to emerge around A = Um- (b) 
For large V extra shot-noise is fed-back into the position variance and the mirror 
thermalizes in a state where the equipartition theorem does not hold, (c) Con- 
tour plot of n vs. K/uJm and "P for A = ujm- Optimal cooling is achieved around 
K ~ 0.2 X LOfn- (d) Fidelity between the mirror and intracavity states in the cooling 
regime as a function of increasing intensity G/uJm with different values of n/ojm = 0.2 
(red line), 0.5 (blue), 1 (green) and 2 (yellow). 

3.3 Readout of the mechanical resonator state 

Eq. flM|) shows that the cavity output is sensitive to the resonator position. 
Therefore, after an appropriate calibration, the cavity output noise power 
spectrum provides a direct measurement of the position noise spectrum Sqiu) 
which, when integrated over a;, yields the value of the position variance (see 
Eq. (I3H])). In many experiments [T6][T7][TM9ppT]l22]l23ll^^ this value is 
employed to estimate the final effective temperature of the cooled resonator 
by assuming energy equipartition {Sp^) ^ {^(f) so that n ~ {^<f) — 1/2. 
However, as we have seen above (see Eqs. P6|) . (1471) . (l53l) . (154|) ). equipartition 
does not generally hold and one should rather estimate {dp^) from ^^(a;) using 
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Eq. (138!) . or directly measure independently the resonator momentum. A dif- 
ferent and more direct way of measuring the resonator temperature, borrowed 
from trapped- ion experiments [88j, has been demonstrated in [25j: if the two 
motional sidebands are well resolved and detected via heterodyne measure- 
ment, the height of the two sideband peaks is proportional to n and to + 1, 
so that one can directly measure the occupancy n from the comparison of the 
two peaks. 

However, one should devise a scheme capable of reconstructing the complete 
quantum state of the resonator from the cavity output light, which is the only 
accessible degree of freedom carrying out information about it. In fact, the full 
reconstruction of the quantum state of the resonator is a necessary condition 
for the unambiguous demonstration of the quantum behavior of the mechanical 
resonator, as for example stationary entanglement, which will be discussed in 
the following. A scheme of this kind has been proposed in [M], based on 
the transfer of the resonator state onto the output field of an additional, fast- 
decaying, "probe" cavity mode. In fact, the annihilation operator of this probe 
cavity mode, Op, obeys an equation analogous to the linearization of Eq. f|T4|) . 

5ap = -{Kp + iAp)6ap + iGpOpSq + y/2K^a'p{t), (57) 

where Kp, Ap, Gp, ap, and a^p{t) are respectively the decay rate, the effective 
detuning, the coupling, the intracavity field amplitude, and the input noise 
of the probe cavity mode. The presence of the probe mode affects the system 
dynamics, but if the driving of the probe mode is much weaker so that Ictpl ^ 
\as\, the back-action of the probe mode on the resonator can be neglected. 
If one chooses parameters so that Ap = Um ^ kp,Gp\ap\, one can rewrite 
Eq. (157|) in the frame rotating at Ap = Um for the slow variables 6d{t) = 
6o(t) exp{iumt} and neglect fast oscillating terms, so to get 

6ap = -Kp6dp + i ^ J' 6b + J2Kpa'p (t), (58) 

where 5b = {i5p + Sq) / a/2 is the mechanical annihilation operator. Finally, if 
Kp ^ Gp\ap\/ \/2, the probe mode adiabatically follows the resonator dynamics 
and one has 

Sap ~ t^^6b + ./^CW- (59) 



K. 



V2 



The input-output relation a""* = ^2Kpbap — a!,"' [79] implies 



C*=^^^5^ + <W, (60) 



showing that, in the chosen parameter regime, the output light of the probe 
mode gives a direct measurement of the resonator dynamics. With an appro- 
priate calibration and applying standard quantum tomographic techniques 
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[89] to this output field, one can therefore reconstruct the quantum state of 

the resonator. 

An alternative way to detect the resonator state by means of state transfer 
onto an optical mode, which does not require an additional probe mode, can 
be devised by appropriately exploiting the strong coupling regime. In this 
second example state transfer is realized in a transient regime soon after the 
preparation of the desired resonator state. One sets the cavity onto resonance 
A = so that the system is always stable, and then strongly increases the 
input power in order to make the coupling G very large, G ^ k, ^loTm- Under 
these conditions, coherent evolution driven by radiation pressure dominates 
and one has state swapping from the mechanical resonator onto the intracavity 
mode in a time tgwap — 7r/2G' so that the cavity mode state reproduces the 
resonator state with a fidelity very close to unity. The fidelity of the swap can 
be computed and reads 



where Vi,V2 are the block matrices in Eq. fl32|) . The resulting fidelity under 
realistic conditions is plotted in Fig. (jSji) as a function of G/uJm for t^/oJm = 
0.2 , 0.5, 1 and 2. One can see that the fidelity is close to unity around the 
optimal cooling regime and that in this regime both the mechanical resonator 
and intracavity field thermalize in the same state. Under this condition one 
can reconstruct the quantum state of the mechanical mode from the detection 
of the cavity output. 



4 Entanglement generation with a single driven cavity mode 

As discussed in the introduction, a cavity coupled to a mechanical degree of 
freedom is capable of producing entanglement between the mechanical and 
the optical modes and also purely optical entanglement between the induced 
motional sidebands. In the following we elucidate the physical origins of this 
entanglement and analyze its magnitude and temperature robustness. More- 
over, we analyze its use as a quantum-communication network resource in 
which the mechanical modes play the role of local nodes that store quantum 
information and optical modes carry this information among nodes. To this 
purpose we apply a multiplexing approach that allows one, by means of spec- 
tral filters, to select many traveling output modes originating from a single 
intracavity field. 



F= y^det (Vi + V2) + (det Vi - 1/4) (det V2 - 1/4) 




(61) 
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4-1 Intracavity optomechanical entanglement 

Entanglement can be easily evaluated and quantified using the logarithmic 
negativity of Eq. (ETi) . which requires the knowledge of the CM of the system of 
interest. For the steady state of the intracavity field- resonator system, the CM 
is determined in a straightforward way by the solution of Eq. f l29|) . However, 
before discussing the general result we try to give an intuitive idea of how 
robust optomechanical entanglement can be generated, by using the sideband 
picture. Using the mechanical annihilation operator 6b introduced in the above 
section, the linearized QLE of Eqs. (fT^ - ffTSl) can be rewritten as 

Sb = -^(6b- 56te2i^-*) + ^b'"' + (5Sl'e^(^+^™)* + 6ae''-'^"^-^\^^) 
6h = -K6a + (55te*(^+-™)* + 6be'^^-^^^^) + v^S*". (63) 

We have introduced the tilded slowly evolving operators 6b{t) = (56(t)e*'^'"*, 
6a{t) = 6a{t)e''^\ and the noises = a*"(t)e^^* and = ^{t)e'''^^/V2. 

The input noise a*"'(t) possesses the same correlation function as a^"'{t), while 
the Brownian noise 6*"()f:) in the limit of large mechanical frequency Um ac- 
quires "optical-like" correlation functions (6*"'"''(t)6*"(t')) = no6{t — t') and 
(6*"(t)6^"'t(t')) = [no + 1] - t') [90J. Eqs. (162D-(163D show that the cavity 
mode and mechanical resonator are coupled by radiation pressure via two 
kinds of interactions: i) a down-conversion process with interaction Hamilto- 
nian 6b'^6a) + 6a6b, which is modulated by a factor oscillating at ujm + A; ii) 
a beam-splitter-like process with interaction Hamiltonian 6b'^6a + 6a^6b, mod- 
ulated by a factor oscillating at — A. Therefore, by tuning the cavity into 
resonance with either the Stokes sideband of the driving laser, A = —ujm, or 
the anti-Stokes sideband of the driving laser, A = Um, one can resonantly 
enhance one of the two processes. In the rotating wave approximation (RWA), 
which is justified in the limit of Um ^ G, k, the off-resonant interaction oscil- 
lates very fast with respect to the timescales of interest and can be neglected. 
Therefore, in the RWA regime, when one chooses A = —uJm, the radiation pres- 
sure induces a down- conversion process, which is known to generate bipartite 
CV entanglement. Instead when one chooses A = uj^, the dominant process 
is the beam-splitter-like interaction, which is not able to generate optome- 
chanical entanglement starting from classical input states [9T], as in this case. 
This argument leads to the conclusion that, in the RWA limit Um ^ G, k, the 
best regime for optomechanical entanglement is when the laser is blue-detuned 
from the cavity resonance A = —uj^ and down-conversion is enhanced. How- 
ever, this argument is valid only in the RWA limit and it is strongly limited by 
the stability conditions, which rather force to work in the opposite regime of a 
red-detuned laser. In fact, the stability condition of Eq. (120|) in the RWA limit 
A = —Um ^ fi^, 7m, simplifies to G < 1/2^7^. Since one needs small mechani- 
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cal dissipation rate 7^ in order to see quantum effects, this means a very low 
maximum value for G. The logarithmic negativity Ej^ is an increasing func- 
tion of the effective optomechanical coupling G (as expected) and therefore 
the stability condition puts a strong upper bound also on It is possible 
to prove that the following bound on Ej\f exists [69] 




showing that < ln2 and above all that entanglement is extremely fragile 
with respect to temperature in the blue-detuned case because, due to the 
stability constraints, Ej^ vanishes as soon as rio > 1- 

This suggests that, due to instability, one can find significant intracavity op- 
tomechanical entanglement, which is also robust against temperature, only far 
from the RWA regime, in the strong coupling regime in the region with positive 
A, because Eq. ( 1221) allows for higher values of the coupling {G < \J + uj'^ 
when A = Um) ■ This is confirmed by Fig. where the exact Ej^ calculated 
from the solution of Eq. fl29l) is plotted versus the normalized detuning A/u;^ 
and the normalized effective optomechanical coupling G/uJm- One sees that 
Ej\f reaches significant values close to the bistability threshold; moreover it is 
possible to see that such intracavity entanglement is robust against thermal 
noise because it survives up to reservoir temperatures around 20 K [51]. It 
is also interesting to compare the conditions for optimal entanglement and 
cooling in this regime where the cavity is resonant with the anti-Stokes side- 
band. In Fig. [6)d, n is plotted versus the same variables in the same parameter 
region. One can see that, while good entanglement is accompanied by good 
cooling, optimal entanglement is achieved for the largest possible coupling 
G allowed by the stability condition. This condition is far from the optimal 
cooling regime, which does not require very large G because otherwise the 
radiation pressure noise contribution and consequently the position variance 
become too large (see Eq. fl54|) and Fig. [5]) [69] . 



4-2 Entanglement with output modes 

Let us now define and evaluate the entanglement of the mechanical resonator 
with the fields at the cavity output, which may represent an essential tool 
for the future integration of micromechanical resonators as quantum memo- 
ries within quantum information networks. The intracavity field 5a{t) and its 
output are related by the usual input-output relation [79j 




Em < In 
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(64) 




(65) 
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Fig. 6. Intracavity entanglement and cooling in the red-detuned regime, (a) Contour 
plot of logarithmic negativity of the field-mirror system at the steady state as a 
function of G/iOrn ^-nd ^j^jJm for the parameters po ^-iid k = uj^i- (b) n in the same 
parameter region: the plot shows that optimal cooling and optimal entanglement 
are both achieved close to A/wm ~ 1. However, optimal cooling is obtained for 
smaller values of G/u)m with respect to entanglement. 

where the output field possesses the same correlation functions of the optical 
input field a*"(t) and the same commutation relation, i.e., the only nonzero 
commutator is a°"*(t), a°"*(t')^ = 5{t — t'). From the continuous output field 
a"^^{t) one can extract many independent optical modes, by selecting different 
time intervals or equivalently, different frequency intervals (see e.g. [92j). One 
can define a generic set of output modes by means of the corresponding 
annihilation operators 



{t)= f dsgk{t-s)a"^\s), k = l, 



(66) 



where gk{s) is the causal filter function defining the A;-th output mode. These 



annihilation operators describe N independent optical modes when 
Sjk, which is verified when 



af*(t),af*(t)1 



dsgj{sygk{s) = 6jk, 



(67) 



i.e., the filter functions gk{t) form an orthonormal set of square-integrable 
functions in [0, oo). The situation can be equivalently described in the fre- 
quency domain: taking the Fourier transform of Eq. fl66l) . one has 
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where gk{^^) is the Fourier transform of the filter function. An exphcit example 
of an orthonormal set of filter functions is given by 

{9 denotes the Heavyside step function) provided that flk and r satisfy the 
condition 

2n 

Qj — fifc = — p, integer p. (70) 

These functions describe a set of independent optical modes, each centered 
around the frequency and with time duration r, i.e., frequency bandwidth 
~ 1/r, since 



When the central frequencies differ by an integer multiple of 27r/r, the cor- 
responding modes are independent due to the destructive interference of the 
oscillating parts of the spectrum. 

The entanglement between the output modes defined above and the mechan- 
ical mode is fully determined by the corresponding (2A^ + 2) x (2A^ + 2) CM, 
which is defined by 

= \ {uT\t)uT\t) + u^^\t)ur\t)) , (72) 

where 

«-*(t) = (5q{t),5p{t),Xl^\t),Yr\t),...,X'^\t),Y°-\t)f (73) 

is the vector formed by the mechanical position and momentum fluctua- 
tions and by the amplitude {X'^\t) = [af*(t) +af*(t)^] 7^2), and phase 

(y^°«*(t) = [af *(t) - af'^ty] Ii\f2) quadratures of the output modes. The 
vector u°'^^{t) properly describes + 1 independent CV bosonic modes, and in 
particular the mechanical resonator is independent of (i.e., it commutes with) 
the N optical output modes because the latter depend upon the output fleld 
at former times only [s < t). From the intracavity CM and Eqs. fl65l) . fl66|) . and 
f l72D one can determine the (A^ + 1) x (A^ + 1) CM matrix V°"* at the steady 
state [69j. 

Let us flrst consider the case when we select and detect only one mode at the 
cavity output. Just to fix the ideas, we choose the mode specified by the filter 
function of Eqs. ( l69l) and ( 17T1) . with central frequency Q and bandwidth r~^. 
Straightforward choices for this output mode are a mode centered either at 
the cavity frequency, f2 = — cu/, or at the driving laser frequency, = (we 
are in the rotating frame and therefore all frequencies are referred to the laser 
frequency a;;), and with a bandwidth of the order of the cavity bandwidth 
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Fig. 7. Resonator-output field entanglement when the central frequency of the 
output mode is swept around the laser frequency. Parameters are po, A = uJm, 
G = uJm/'^ and k = uj^n- The entanglement is optimized when the output mode 
coincides with the Stokes sideband of the laser ($7 = —uJm), with the appropriate 
bandwidth (e ~ 10, corresponding to Tjm'^ ~ !)• For smaller e, the selected output 
mode mixes Stokes and anti-Stokes photons and the entanglement is weak, while 
for larger e only a fraction of the sideband is selected and part of the quantum 
correlations are lost. In the inset the robustness of Stokes-mirror Ej\/ with respect 
to temperature is shown. 

r^^ K. However, as discussed above, the motion of the mechanical resonator 
generates Stokes and anti-Stokes motional sidebands, consequently modifying 
the cavity output spectrum. 

In order to determine the output optical mode which is better entangled with 
the mechanical resonator, we study the logarithmic negativity Ej\/ associated 
with the output CM V°"* (for = 1) as a function of the central frequency of 
the mode fl and its bandwidth r~^, at the same parameter region considered 
in the previous subsection, pq and A = Um, where intracavity entanglement 
is optimal. The results are shown in Fig. [7l where Ej^ is plotted versus fl/um 
at different values of e = rUm- If ^ 1, i-e., the bandwidth of the detected 
mode is larger than Um, the detector does not resolve the motional sidebands, 
and Ef^f has a value (roughly equal to that of the intracavity case) which 
does not essentially depend upon the central frequency. For smaller band- 
widths (larger e), the sidebands are resolved by the detection and the role of 
the central frequency becomes important. In particular Ej^ becomes highly 
peaked around the Stokes sideband Q = —utm, showing that the optomechan- 
ical entanglement generated within the cavity is mostly carried by this lower 
frequency sideband. What is relevant is that the optomechanical entanglement 
of the output mode is significantly larger than its intracavity counterpart and 
achieves its maximum value at the optimal value e ~ 10, i.e., a detection band- 
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width ^ uJm/10. This means that in practice, by appropriately filtering 
the output light, one realizes an effective entanglement distillation because the 
selected output mode is more entangled with the mechanical resonator than 
the intracavity field. 

The fact that the output mode which is most entangled with the mechanical 
resonator is the one centered around the Stokes sideband is also consistent 
with the physics of a previous model analyzed in [BBj. In [66j, a free-space 
optomechanical model is discussed, where the entanglement between a vibra- 
tional mode of a perfectly refiecting micro-mirror and the two first motional 
sidebands of an intense laser beam shined on the mirror is analyzed. Also in 
that case, the mechanical mode is entangled only with the Stokes mode and 
it is not entangled with the anti-Stokes sideband. 

One can also understand why the output mode optimally entangled with the 
mechanical mode has a finite bandwidth ^ tUm/lO (for the chosen op- 
erating point). In fact, the optimal situation is achieved when the detected 
output mode overlaps as best as possible with the Stokes peak in the spec- 
trum, and therefore coincides with the width of the Stokes peak. This 
width is determined by the effective damping rate of the mechanical resonator, 
1m ^ = 7m + r, given by the sum of the intrinsic damping rate 7^ and the net 
laser cooling rate F of Eq. fl52p . It is possible to check that, with the chosen 
parameter values, the condition e = 10 corresponds to r^^ ~ 7^-^. 

It is finally important to analyze the robustness of the present optomechani- 
cal entanglement with respect to temperature. As discussed above and shown 
in [M], the entanglement of the resonator with the intracavity mode is very 
robust. It is important to see if this robustness is kept also by the optome- 
chanical entanglement of the output mode. This is shown also in the inset of 
Fig. [71 where the logarithmic negativity Ej^ of the output mode centered at 
the Stokes sideband Q = —oj^ is plotted versus the temperature of the reser- 
voir at two different values of the bandwidth, the optimal one e = 10, and at 
a larger bandwidth e = 0.5. We see the expected decay of Ej^ for increasing 
temperature, but above all that also this output optomechanical entangle- 
ment is robust against temperature because it persists even above liquid He 
temperatures, at least in the case of the optimal detection bandwidth e = 10. 



4-3 Optical entanglement between sidebands 

Let us now consider the case where we detect at the output two independent, 
well resolved, optical output modes. We use again the step-like filter functions 
of Eqs. (l69l) and (!7T!) . assuming the same bandwidth for both modes 
and two different central frequencies, fli and fl2, satisfying the orthogonality 
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Fig. 8. Sideband-sideband entanglement. Parameters po, k = oJrn and G = a'm/2. 
(a) Assuming one detection setup centered at the Stokes sideband and sweeping 
the second detection frequency around the anti-Stokes sideband at O = ujm, the 
entanglement is clearly shown to be optimized when the anti-Stokes output field 
is detected. This entanglement is improving with smaller and smaller detection 
bandwidth (e — > oo). (b) Logarithmic negativity increases asymptotically to a finite 
value with e oo. (c) Temperature robustness for e = IOtt and e = IOOtt. The 
entanglement survives to very high temperatures. 

condition of Eq. ( |70l) VL1 — VL2 = 2p7rr~^ for some integer p, in order to have two 
independent optical modes. It is interesting to analyze the stationary state of 
the resulting tripartite CV system formed by the two output modes and the 
mechanical mode, in order to see if and when it is able to show purely optical 
bipartite entanglement between the two output modes. 



The generation of two entangled light beams by means of the radiation pres- 
sure interaction of these fields with a mechanical element has been already con- 
sidered in various configurations. In Ref. [93], and more recently in Ref. |63j, 
two modes of a Fabry-Perot cavity system with a movable mirror, each driven 
by an intense laser, are entangled at the output due to their common pondero- 
motive interaction with the movable mirror (the scheme has been then gen- 
eralized to many driven modes in p4]). In the single mirror free-space model 
of Ref. [66], the two first motional sidebands are also robustly entangled by 
the radiation pressure interaction as in a two-mode squeezed state produced 
by a non-degenerate parametric amplifier ^95j. Robust two-mode squeezing of 
a bimodal cavity system can be similarly produced if the movable mirror is 
replaced by a single ion trapped within the cavity [96] . 

The situation considered here is significantly different from that of Refs. [63|93]|9^ 
which require many driven cavity modes, each associated with the correspond- 
ing output mode. In the present case instead, the different output modes orig- 
inate from the same single driven cavity mode, and therefore it is simpler 
from an experimental point of view. The present scheme can be considered 
as a sort of "cavity version" of the free-space case of Ref. [66], where the 
reflecting mirror is driven by a single intense laser. Therefore, as in [66|95] . 
one expects to find a parameter region where the two output modes centered 
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around the two motional sidebands of the laser are entangled. This expec- 
tation is clearly confirmed by Fig. [8^, where the logarithmic negativity Ej^ 
associated with the bipartite system formed by the output mode centered at 
the Stokes sideband {Vti = —Um) and a second output mode with the same 
inverse bandwidth {e = uJmT = lOvr) and a variable central frequency fi, is 
plotted versus VL/uJm- Ej^f is calculated from the CM V°"* (for N = 2) eliminat- 
ing the first two rows associated with the mechanical mode. One can clearly 
see that bipartite entanglement between the two cavity outputs exists only 
in a narrow frequency interval around the anti-Stokes sideband, Q = Um, 
where Ej^ achieves its maximum. This shows that, as in [SSES], the two cav- 
ity output modes corresponding to the Stokes and anti-Stokes sidebands of 
the driving laser are significantly entangled by their common interaction with 
the mechanical resonator. The advantage of the present cavity scheme with 
respect to the free-space case of p^ll^ is that the parameter regime for reach- 
ing radiation-pressure mediated optical entanglement is much more promising 
from an experimental point of view because it requires less input power and 
a not too large mechanical quality factor of the resonator. In Fig. [8)d, the 
dependence of Ej^ of the two output modes centered at the two sidebands 
f2 = ±uJm upon their inverse bandwidth e is studied. We see that, differently 
from optomechanical entanglement of the former subsection, the logarithmic 
negativity of the two sidebands always increases for decreasing bandwidth, 
and it achieves a significant value, comparable to that achievable with para- 
metric oscillators, for very narrow bandwidths. This fact can be understood 
from the fact that quantum correlations between the two sidebands are estab- 
lished by the coherent scattering of the cavity photons by the oscillator, and 
that the quantum coherence between the two scattering processes is maximal 
for output photons with frequencies ui ± Um- Figs. [7] and [8] show that in the 
chosen parameter regime, the output mode centered around the Stokes side- 
band mode shows bipartite entanglement simultaneously with the mechanical 
mode and with the anti-Stokes sideband mode. This fact suggests that the CV 
tripartite system formed by the output Stokes and anti-Stokes sidebands and 
the mechanical resonator mode might be characterized by a fully tripartite- 
entangled stationary state. This is actually true and it can be checked by 
applying the classification criterion of Ref . [97] , providing a necessary and suf- 
ficient criterion for the determination of the entanglement class in the case of 
tripartite CV Gaussian states, which is directly computable in terms of the 
eigenvalues of appropriate test matrices [97] (see Ref. [69]). 



5 Entanglement generation with two driven cavity modes 

We now generalize the system by considering the case when two cavity modes 
with different frequencies are intensely driven. We shall focus onto a param- 
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eter regime which will prove to be convenient for the generation of robust 
stationary CV entanglement between the resonator and the two cavity modes. 
A bichromatic driving of a cavity has been already experimentally considered 
in Refs. [2I|. There however it was employed for cooling a macroscopically 
heavy (m ~ Ig) movable mirror. One driven mode is used to "trap" the mir- 
ror, i.e., to induce a strong optical spring effect increasing by three orders of 
magnitude the oscillation frequency. The other driven mode is instead used to 
cool the mechanical resonator by increasing the effective mechanical damping, 
either via back-action, or via cold-damping feedback. The bichromatic driving 
configuration has been already considered for the generation of entanglement 
in various configurations in some theoretical proposals. In fact, in Ref. |93j . 
and more recently in Ref. [63], two modes of a Fabry- Perot cavity system, each 
driven by an intense laser, are entangled at the output due to their common 
ponderomotive interaction with the movable mirror. 

5.1 Quantum Langevin equations and stability conditions 

We generalize the Hamiltonian of Eq. by considering two cavity modes, 
with frequencies Uca and uJcB, each driven by a laser with frequency Uqa and 
ivqb, and power Va and Vb, respectively. The resulting Hamiltonian is 

H = hucA a^CL + huJcB b^b + -hujrn{p^ + q^) — ^{Gqa Q^Of + Gqb b^b)q 
+ih[EA{ah-'''''^' - ae^'^o-^*) + Esib^-''^^^' - 6e*"°^*)], (74) 

where a and b now denote the annihilation operators of the two cavity modes, 
we have introduced the two coupling constants Gq^ = ^ h/mum uJcx/L, and 

the two driving rates \Ex\ = \j2PxK,/tujjQxi x = A,B. We have assumed for 
simplicity that the two modes have the same decay rate k. We assume that 
scattering of photons of the driven modes into other cavity modes and also 
between the two chosen modes is negligible: this is verified when Um is much 
smaller then the free spectral range of the cavity. 

Introducing again dissipation and noise terms as in Sec. [2], the system dynamics 
is described by the following set of nonlinear QLE, written in the interaction 
picture with respect to /zo^oao^o, + huoBb^b, 



q = UmP, (75) 

P = -Umq - ImP + GqaO^O + G'oB&"^& + ^, (76) 

a = -[K + i(AoA - GQAq)]a + Ea + v/2^a*", (77) 

b = -[K + ^(AoB - GoBq)]b + Eb + V2^b'^, (78) 
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where Aqx = ^cx — t^ox are the detunings of the two lasers, and we have 
introduced a vacuum input noise for the cavity mode 6, possessing the 

same correlations of Eqs. ([7j)-([HD- 

We assume again that both modes are intensely driven so that the system is 
characterized by a semiclassical steady state with large intracavity amplitudes 
for both modes and a modified cavity length. This classical steady state is 
determined by setting the time derivatives to zero, factorizing the averages 
and solving for the mean values = (a), = (&), Qs = {q), Ps = (p)- One 
gets 



Ea 



Qs 

Ps 



where the effective detunings = Aq^ — {G\ 



0A\ 



X 



(79) 
(80) 

(81) 
(82) 



have been defined, so that Eqs. fl79l) - fl80l) form actually a system of nonlinear 
equations, whose solution gives the stationary amplitudes and bg. 

One then focuses onto the dynamics of the quantum fluctuations around this 
steady state, which are well described by linearizing the QLE of Eqs. f l75l) - fl78l) 
around the semiclassical steady state values, provided that \as\, \bs\ ^ 1. The 
linearized QLE for the resonator and for the amplitude and phase quadratures 
of the two modes, 6Xa, SXb, SYa and 6Yb, defined as in Sec. Ej can be written 
in compact form as 

u{t) = Au{t) + n{t), 

where u = {6q, 6p, SXa, SYa, SXb, SYb)^ is the vector of quadrature fluctua- 
tions, and n = (0,^, V^Xj', V2kYX\ V^kX'^, ^/2kY}^Y is the correspond- 
ing vector of noises. The 6x6 matrix A is the drift matrix of the system, 
which reads 

^Ocj^O 0^ 
-u^ 7„ Gb 

-K Aa 
Ga -Aa 

-K Ab 
\ Gb -Ab ) 
where we have chosen the phase reference of the two cavity modes so that a 



A 



83) 
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and bg are real and positive, and defined the effective couplings Ga = Gqa^sV^ 
and Gb = GoBbsV^- 

The steady state exists and it is stable if all the eigenvalues of the drift ma- 
trix A have negative real parts. The parameter region under which stability 
is verified can be obtained from the Routh-Hurwitz criteria [81], but the in- 
equalities that come out in the present case are quite involved. One can have 
an idea of this fact from the expression of the characteristic polynomial of A, 
P(A) = + ciA^ + CsA^ + CgA^ + C4A2 + C5A + C6, where 



ci = 7m + 4k, 

C2 = + A| + A-fmK + Gk'^ + ul, 

C3 = 7m(AA + A| + + 2k[AI + A| + + fi^)], 

C4 = K^ + 27^/€(A| + 2k^) + Qk^ujI, + AI{k^ + cu^) + 

A^(A| + 27„/t + K^+ UjI^) - UJra{G\AA + G\Ab), 
C5 = 7m(A^ + /€')(A| + h") + 2/€C^^(A^ + A^ + 21^) 

-2kuj^{G\Aa\ G\Ab), 
C6 = u:1{A\ + «:2)(A| + h") - ujJ,G\Ab{^\ + 

-rG\AA{A\^^)\. 

We are considering here a bichromatic driving of the cavity in order to im- 
prove the size and the robustness of the generated entanglement. Entangle- 
ment monotonically increases with the optomechanical coupling but, as we 
have seen also in the previous sections, the stability conditions put a strict 
upper bound on the maximum achievable value of this coupling. Therefore 
it is interesting to find a regime in which the presence of the second driven 
mode makes the system always stable, so that the couplings can be made 
very large (for example by increasing the input power, the cavity finesse, or 
decreasing the cavity length) without entering the unstable regime. One then 
hopes that in this regime also entanglement can be made large and robust 
against temperature. 

A simple way to have always a stable system is to find a particular relation 
between the parameters such that the characteristic polynomial of A does 
not depend upon Ga and Gb- In this case, the eigenvalues of A would be 
independent of the two couplings and stability would be guaranteed. The 
expressions above show that the eigenvalues of A are independent of Ga and 
Gb and the system is always stable if and only if 



\Ga\ = \Gb\ = G, (84a) 
Aa = -Ab = A. (84b) 

The condition described by Eqs. represents a perfect balance between a 
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cooling cavity mode (which, without loss of generality, we can take as mode 
A, so that > 0) and a heating cavity mode, i.e., mode B with < 0. The 
fact that the eigenvalues of A do not depend upon the couplings means that 
the decay rates of both the resonator and the cavity modes are left unchanged 
and in this case radiation pressure mainly create quantum correlations, i.e., 
entanglement, between the modes. We shall assume conditions flH^ from now 
on. 



5.2 Entanglement of the output modes 

We now calculate the entanglement properties of the steady state of the bichro- 
matically driven cavity. However we shall not discuss here the intracavity en- 
tanglement, but only the entanglement properties of the optical output modes. 
In fact, as we have seen above in the case of a single driven mode, one can ob- 
tain a larger optomechanical entanglement with respect to the intracavity case 
by appropriately filtering the output modes. Moreover only the entanglement 
with output modes is relevant for any quantum communication application. 
We shall apply therefore the filter function formalism developed in Sec. 14. 2[ 
restricted however here to the simple case of a single output mode for each 
intracavity mode. In fact, we have now two driven cavity modes and consider- 
ing the more general case of multiple output modes associated to each driven 
mode as in Sec. 14.21 would render the description much more involved with- 
out however gaining too much insight into the physics of the problem. The 
two output modes originate from two different cavity modes, and since the 
latter are not too close in frequency, they consequently describe two indepen- 
dent modes. Therefore we do not need orthogonal filter functions like those of 
Eq. (1691) used for the single driven mode case, and we choose here a different 
filter function. We consider the two output modes with annihilation operators 



where a°^^{t) and a°j!^^{t) are the usual output fields associated with the two 
cavity modes and 



are the two filter functions, describing two output modes, both with bandwidth 
1/r and with central frequencies, VLa and VLb., which are in general different 
from the cavity mode frequencies oJca and uJcb- 

The entanglement between the chosen output modes and the mechanical res- 
onator mode is fully determined by the corresponding 6x6 CM, which is 




(85) 




(86) 
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defined as in Eq. (172|) 



where now 



w°"*(t) = 0, 0, 5x-;(t), 5X-*(t), 6Y^-Ht) 



is the vector formed by the mechanical position and momentum fluctuations 
and by the amplitude and phase quadratures of the filtered modes. Using the 
various definitions, input-output relations and also the correlation function of 
the noise terms, one can derive an integral expression for the CM V°"* of the 
system (see Ref. [21] for the details in a similar calculation), which is given by 
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where T{uj) is the Fourier transform of 
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(90) 



M{uj) = {iuj + Ay\ (91) 

Pout = Diag[0,0, 1, 1, 1, 1] is the projector onto the optical quadratures, and 
D{uj) is the matrix associated with the Fourier transform of the noise corre- 
lation functions, given by 

D{uj) = Diag[0, {'-fm^ / ujjn) cothlhw/ 2k bT), 

Using the CM V°"* one can analyze the entanglement between the three dif- 
ferent bipartitions of the system, when one of the three modes is traced out, 
and also tripartite entanglement. 



5.2.1 Optomechanical entanglement 

First of all we consider the entanglement between the output field of the 
"cooling mode" (A) (the one with > 0) and the mechanical resonator. 
We have seen in Sec. IHthat this configuration allows to achieve the maximum 
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optomechanical entanglement in the case of a single driven cavity mode. In 
fact, when A ~ tUm, G is large enough, and the selected output mode is 
centered around the Stokes sideband, the entanglement is optimized and it is 
also robust against temperature (see Fig. [7]). Fig. [3 shows that the presence 
of the second "heating" mode B disturbs this optimal condition and that Ej^ 
is appreciably lower than the one with only one driven mode. In fact, we 
have considered here a similar parameter region, i.e. po, ^^ = f^m, = (^m, 
Ab = —uJm, Ga = 0.326a;m, = 0.302u;m. The qualitative behavior of Ej^ 
is identical to that of the corresponding Fig. [TJ i.e., Ej^ is maximum when 
the output mode overlaps as best as possible with the Stokes sideband of the 
corresponding driving laser, which means centered around —Um and with an 
inverse bandwidth e = ~ 10. However the achievable values of E_\f are 
significantly lower. Fig. [7)d shows that, despite the lower values, entanglement 
is still quite robust against temperature. 




Fig. 9. Logarithmic negativity Ej\f of the bipartite system formed by the mechanical 
mode and the output of the "cooling" mode A. (a) Ej\f versus the center frequency of 
the output mode O^/wm at three different values of the inverse detection bandwidth 
e = LOmT- As in the single driven mode case (see Fig. [7]), entanglement is maximum 
when the output mode is centered around the the Stokes sideband ^Ia_ = —^m- 
The other parameters are pq, n = ujm-, = ^m-, = — Ga = 0.326a;m) 
Gh = 0.302u;m- (b) Ej^f versus the reservoir temperature T when the output mode 
is centered at the Stokes sideband {VIa = —^m) for the same three different values 
of e. 

The advantage of the bichromatic driving becomes instead apparent when 
one considers the bipartite system formed by the resonator and the output 
field of the "heating" mode (B), the one with = —^j^m- The stationary 
optomechanical entanglement can achieve in this case significantly larger val- 
ues. The results are shown in Fig. [10] which refers to the same parameters of 
Fig. [9] and shows the same qualitative behavior: Ej^ is optimized when the 
selected output mode well overlaps with the Stokes sideband of the driving 
laser VLb = —oJm and it persists up to reservoir temperatures of the order of 
10 K. However, Ej^ is now roughly three times larger than the corresponding 
value for the "cooling" mode. This behavior is different from what is found in 
Sec. m for a single driven cavity mode, where we have seen that optomechan- 
ical entanglement in the "heating" regime of negative detunings is seriously 
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Fig. 10. Logarithmic negativity Ej\/ of the bipartite system formed by the mechanical 
mode and the output of the "heating" mode B. (a) Ej^ versus the center frequency 
of the latter ^a/^tti at three different values of the inverse detection bandwidth 
e = LOmT- As it happens for the "cooling" mode A, entanglement is maximum when 
the output mode is centered around the the Stokes sideband = —^m- Parameters 
are as in Fig. [9l (b) Ej\f versus the reservoir temperature T when the output mode 
is centered at the Stokes sideband (f^B = —ujm) for the same three different values 
of e. 

limited by stability conditions. Now, thanks to the combined action of the two 
driven modes and to the conditions fl84l) . the system is always stable and the 
parametric-like process described in Sec. H] is able to generate large and robust 
entanglement. Therefore we can say that in this bichromatic case, mode A 
helps to entangle in a robust way the output of the "heating" mode B, by 
counteracting its instability effect and making the system stable for any value 
of the couplings Ga and Gb- Notice that in this case, the Stokes sideband 
of the laser driving mode B is resonant with the cavity, because = —Um 
implies oJcb = ^ob — = ^stokes and this provides a further reason why the 
optomechanical entanglement may become large. 

5.2.2 Purely optical entanglement between output modes 

Let us consider now the purely optical entanglement between the two output 
hght beams. As discussed at the beginning of the section, the possibility to 
entangle two different output modes of a cavity by means of radiation pres- 
sure has been already suggested in different configurations [63]I66II93|95] . We 
have also seen in Sec. H] that this is possible even with a single driven mode. 
It is nonetheless interesting to compare the results of Sec. Hlwith the present 
bichromatic driving case. The bichromatic case has been already studied in 
ReL [63], which however restricted to the case of output modes with infinitely 
narrow bandwidth (r = oo) and centered around the driving laser frequency 
{VLa = fi^ = 0). The general filter function formalism instead allows us to con- 
sider generic values of r, and fi^. By applying again Eq. (189|) and tracing 
out now the mechanical mode, we get the results illustrated in Fig. [TTl We have 
considered a slightly different parameter regime with respect to the previous 
subsection, by choosing slightly larger couplings, Ga = 1.74a;m, Gi, = l.TOwm, 
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Fig. 11. Logarithmic negativity Ej^ of the bipartite system formed by the output 
modes associated with the two driven cavity modes, (a) Ej^ versus the center fre- 
quency of the "coohng" mode A Qjn/um for a center frequency of the "heating" mode 
fixed at ^Ib = —^m (Stokes sideband), and at three different values of the inverse 
detection bandwidth e = UmT- The other parameters are po, k = ojm-, = ^m-, 
Ab = —oJmi Ga = l-74u;m> G}, = l.TOwm- (b) Ej^f versus the reservoir temperature T 
when the output of the mode A is centered at the anti-Stokes sideband {Q,a = oJm) 
and the output of mode B is centered at the Stokes sideband {^b = —^m), for two 
different values of e. 

i.e., larger input powers. Here, the oscillating mirror induces Stokes and anti- 
Stokes sidebands for both driving lasers and therefore it may be nontrivial to 
establish which are the most-entangled output modes. Fig. [TlT a) shows that 
the largest all-optical entanglement is achieved between the anti-Stokes side- 
band of the "cooling" mode and the Stokes sideband of the "heating" beam. 
This is consistent with the results for a single cavity mode, because in both 
cases the motion of the resonator creates strong quantum correlations between 
the scattering of a Stokes and an anti-Stokes photon. Moreover this result can 
be understood from the fact that the two sidebands are those which are res- 
onant with the corresponding cavity mode. Fig. [TlT a) also shows that, as it 
happens in the single cavity mode case, and differently form the optomechan- 
ical entanglement, the all-optical E/^ monotonically increases for decreasing 
detection bandwidths. This is reasonable because the two output modes are 
correlated as in two-mode squeezing which is based on the pairwise correlated 
production of photons from a pump laser beam via a parametric process. In 
this case the quantum correlations are optimally detected when only pairs of 
photons exactly satisfying the matching condition Ug + ujas = ^oa + ^^ob are 
detected, i.e., when t = oo. 

Fig. fTlT b) instead shows the robustness of all-optical entanglement with re- 
spect to the reservoir temperature, which is extremely good: entanglement 
persists even at room temperature provided that one considers output modes 
with a sufficiently narrow bandwidth. In this respect, the bichromatic driving 
case proves to be more promising than the single driving mode case (compare 
Fig.imb) with Fig.E^c)). 
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Combining all the results of this section, we see that the output modes as- 
sociated with the two driven cavity modes and the mechanical mode form a 
tripartite system in which each bipartite subsystem is entangled. This suggests 
that a parameter region exists where this tripartite system is characterized by 
a fully tripartite-entangled stationary state. This is actually true and it can be 
checked by applying the classification criterion of Ref. [9^, providing a nec- 
essary and sufficient criterion for the determination of the entanglement class 
in the case of tripartite CV Gaussian states, which is directly computable in 
terms of the eigenvalues of appropriate test matrices ^7\. 



6 Cavity-mediated atom-mirror stationary entanglement 

A final recent application of optomechanical systems, recently suggested by a 
number of papers (see Refs. [70H7T] ). is to couple them also to atomic ensembles 
in order to realize new and more fiexible CV quantum interfaces. To be more 
specific, here we consider a hybrid system comprised of A^^ two-level atoms 
of energy splitting hwa, coupled to an optical cavity, which is in turn coupled 
to a mechanical element by radiation pressure. We consider again the steady 
state of the system and choose a weak-coupling regime where the atoms and 
the cavity are far-off resonance (as illustrated by Fig. [T2|) . The working point 
for the optomechanical system is the regime described in the previous section 
where red-detuned driving of the cavity ensures optimal entanglement between 
the Stokes sideband and the mechanical resonator. We show here that when 
the atoms are resonant with the Stokes sideband of the laser, a regime where 
both atoms-mirror bipartite CV entanglement and tripartite CV entanglement 
can be generated in the steady state, is achieved. 

We start from the Hamiltonian of Eq. ([T]) to which we add the Tavis-Cummings 
atom-cavity field interaction 

Hj = hg (^S+a + S-a^^ , 

where collective spin operators are defined as S+^^^z = o"+!-,z for ^ = 
l,Na {a+^-^z are the Pauli matrices) and satisfy the commutation relations 
[5*+, S*.] = 5*2 and [5*2, S-t] = ±25*31-. The atom-cavity coupling constant is 

given hj g = uJc/'^heoV where V is the cavity mode volume, is the dipole 
moment of the atomic transition, and eo is the free space permittivity. 

The djTiamics of the tripartite system is fairly complicated. However, one can 
find a regime where a simpler dynamics of three coupled harmonics oscillators 
is a good approximation of the system dynamics. To this purpose, we assume 
that the atoms are initially prepared in their ground state, so that Sz — 
(Sz) — —Na and this condition is not appreciably altered by the interaction 
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Fig. 12. Setup for tripartite hybrid entanglement. An atomic cloud of two- level 
atoms is placed inside a cavity driven by a laser. As seen in the inset, the atoms 
are resonant with the Stokes sideband of the laser. Since this latter sideband is the 
one carrying most of the optomechanical entanglement, also the atoms and movable 
mirror become entangled at the steady state. 

with the cavity field. This is satisfied when the excitation probability of a 
single atom is small. In this limit the dynamics of the atomic polarization can 
be described in terms of bosonic operators: in fact if one defines the atomic 



annihilation operator c = S-/ ^\{Sz)\, one can see that it satisfies the usual 
bosonic commutation relation [c, c^] = 1 [98]. In the frame rotating at the laser 
frequency uji for the atom-cavity system, the quantum Langevin equations can 
then be written as 



where Aq = Uc — oji and = — uji are the cavity and atomic detuning with 
respect to the laser, Ga = g\^Na, and 2'ja is the decay rate of the atomic excited 
level. The noise affecting the atoms has one non-vanishing correlation function 
{Fc {t) {t')) = 6 {t — t') . We now assume that the cavity is intensely driven, 
so that at the steady state, the intracavity field has a large amplitude a^, with 
\as\ S> 1. However, the single-atom excitation probability is g'^\as\'^ / {A'l + •y^) 
and since this probability has to be much smaller than one for the validity 
of the bosonic description of the atomic polarization, this imposes an upper 
bound to Therefore the two conditions are simultaneously satisfied only 
if the atoms are weakly coupled to the cavity, g^/[Al + 7^] <^ Ic^sl^^ ^ 1- 

If one is interested only in atoms- mirror entanglement, one could assume a bad 
cavity limit and adiabatically eliminate the cavity mode [7T]. However, one 




q = UmP, 

p = -uJmQ - ImP + Goo^a + ^, 
a = — (k + 2Ao)a + iGoaq — iGaC 




+ Si + V2Kai, 



(92) 
(93) 
(94) 
(95) 
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can have a more complete information by linearizing the Langevin equations 
Eqs. (!92!)- (p5l) around the semiclassical steady state and then solving for the 
exact solution of the 3-mode system steady state provided by the Lyapunov 
equation fl29p [70j . In fact, owing to the Gaussian nature of the quantum noise 
terms ^, Qin and Fc, and to the linearization of the dynamics, the steady state 
of the quantum fluctuations of the system is a CV tripartite Gaussian state, 
which is completely determined by its 6 x 6 correlation matrix (CM). The 
corresponding drift matrix after linearization is given by 
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(96) 



while the diffusion matrix is equal to D =diag[0, 7m (2no + !),/€, k, 7^, 7a 1. We 
have solved Eq. fl29|) for the CM V in a wide range of the parameters G, Ga, A 
and Aq. We have studied first of all the stationary entanglement of the three 
possible bipartite subsystems, by quantifying it in terms of the logarithmic 
negativity of bimodal Gaussian states. We will denote the logarithmic neg- 
ativities for the mirror-atom, atom-field and mirror-field bimodal partitions 
with Ema, Eaf and Emf, respectively. 

The results on the behavior of the bipartite entanglement Ema are shown in 
Fig. [TSk . Optimization requires, as expected that the atoms are resonant with 
the Stokes motional sideband. In Fig. [JSt, the logarithmic negativity of the 
three bipartitions is plotted versus the normalized atomic detuning. It is evi- 
dent that one has a sort of entanglement sharing: due to the presence of the 
atoms, the initial cavity-mirror entanglement (represented by the dashed line) 
is partially redistributed to the atom-mirror and atom-cavity subsystems and 
this effect is predominant when the atoms are resonant with the Stokes side- 
band {Aa = —uJm)- It is remarkable that, in the chosen parameter regime, the 
largest stationary entanglement is the one between atoms and mirror which are 
only indirectly coupled. Moreover, the nonzero atom-cavity entanglement ap- 
pears only thanks to the effect of the mirror dynamics because in the bosonic 
approximation we are considering and with a fixed mirror, there would be 
no direct atom-cavity entanglement. We also notice that atom- mirror entan- 
glement is instead not present at A^ = Um- This is due to the fact that the 
cavity- mirror entanglement is mostly carried by the Stokes sideband and that, 
when Aa = ujm, mirror cavity-cooling is disturbed by the anti-Stokes photons 
being recycled in the cavity by the absorbing atoms. 
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Fig. 13. Entanglement in the hybrid mirror-atom-field system. Parameters are po, 
K = 7(1 = ujm, G = l.SuJm- (&) Contour plot of Ej\f between mirror and atoms as a 
function of Ga/oOm and Aa/um- The entanglement is optimized for = —oom, i-e. 
when the atoms are resonant with the Stokes sideband of the laser, (b) The three 
bipartite entanglement versus the atomic detuning. The blue dashed line represents 
the mirror-field Ej\f in the absence of atom-field coupling. When the atoms are 
immersed in the mirror-field system, the entanglement is redistributed among the 
three sub-partitions, especially around the regime where Aq = — Wm- 

We notice that the chosen parameters correspond to a small cavity mode vol- 
ume {V ~ 10~^^ m'^), implying that for a dipole transition, g is not small. 
Therefore the assumed weak coupling condition (/^/[A^ + 7a] l^^sl^^ ^ 1 
can be satisfied only if g represents a much smaller, time averaged, coupling 
constant. This holds for example for an atomic vapor cell much larger than 
the cavity mode: if the (hot) atoms move in a cylindrical cell with axis or- 
thogonal to the cavity axis, with diameter ~ 0.5 mm and height ~ 1 cm, they 
will roughly spend only one thousandth of their time within the cavity mode 
region. This yields an effective g ~ 10*^ Hz, so that the assumptions made 
here hold, and the chosen value Ga/'2^T^ = 6 x 10^ Hz can be obtained with 
Na ~ 10^. An alternative solution could be choosing a cold atomic ensemble 
and a dipole-forbidden transition. 

The entanglement properties of the steady state of the tripartite system can be 
verified by experimentally measuring the corresponding CM. This can be done 
by combining existing experimental techniques. The cavity field quadratures 
can be measured directly by homodyning the cavity output, while the mechan- 
ical position and momentum can be measured with the schemes discussed in 
Sec. 13.31 Finally, the atomic polarization quadratures x and y (proportional 
to Sx and Sy) can be measured by adopting the same scheme of Ref. [99], i.e., 
by making a Stokes parameter measurement of a laser beam, shined transver- 
sal to the cavity and to the cell and off-resonantly tuned to another atomic 
transition. 
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7 Conclusions 



The search for experimental demonstrations of the quantum behavior of macro- 
scopic mechanical resonators is a fastly growing field of investigation. Re- 
cent experimental results [TCTlM71[TM9ll20ll2T]|22]^^ 

and theoretical predictions suggest that quantum states of resonators with a 
mass at the microgram scale will be generated and detected in the near future. 
In this chapter we have tried to give an overview of two relevant arguments of 
this research field: i) cooling to the motional ground state; ii) the generation 
of robust entangled steady states involving mechanical and optical degrees 
of freedom. The latter condition is the fundamental prerequisite for the even- 
tual integration of micro- and nano-mechanical resonators serving as quantum 
memories and interfaces within quantum communication networks. 

In the first part of the chapter we have described and compared the two 
main approaches for cooling micro-mechanical resonators via radiation pres- 
sure coupling to an optical cavity, cold-damping feedback [32ll34|35|36|41j , and 
back-action cooling |33|39ll4U|41|42|43j . We have adopted a general quantum 
Langevin treatment which is valid within the full parameter range of a sta- 
ble cavity. Both back-action cooling and cold damping feedback are able to 
cool to the ground state, even though back-action cooling is preferable for a 
good cavity (k < Um) , while cold damping is more convenient for a bad cavity 

In the second part of the chapter we have analyzed the entanglement proper- 
ties of the steady state of the system formed by the optical cavity coupled to 
a mechanical element. We have considered two different configurations, with 
either one or two intensely driven cavity modes. We have seen that the intra- 
cavity mode and the mechanical element can be entangled in a robust way 
against temperature, and that back-action cooling is not a necessary condi- 
tion for achieving entanglement. In fact, entanglement is possible also in the 
opposite regime of a blue-detuned laser where the cavity mode drives and 
does not cool the resonator. More generally, the two phenomena are quite in- 
dependent, and one is not necessarily accompanied by the other. Cooling is a 
classical process (even though it can ultimately lead to the quantum ground 
state), while entanglement is an intrinsically quantum phenomenon. More- 
over, they are optimized in different parameter regimes. In fact, logarithmic 
negativity is maximized close to the stability threshold of the system, where 
instead the resonator is not cooled. We have then focused our study onto the 
entanglement properties of the cavity output field, which is the relevant one 
for quantum communication applications. We have developed a general theory 
showing how it is possible to define and evaluate the entanglement properties 
of the multipartite system formed by the mechanical resonator and N inde- 
pendent output modes of the cavity field. We have seen that the tripartite 
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system formed by the mechanical element and the two output modes centered 
at the first Stokes and anti-Stokes sideband of the driving laser (where the 
cavity output noise spectrum is concentrated) shows robust fully tripartite 
entanglement. In particular, the Stokes output mode is strongly entangled 
with the mechanical mode and shows a sort of entanglement distillation be- 
cause its logarithmic negativity is significantly larger than the intracavity one 
when its bandwidth is appropriately chosen. In the same parameter regime, 
the Stokes and anti-Stokes sideband modes are robustly entangled, and the 
achievable entanglement in the limit of a very narrow detection bandwidth is 
comparable to that generated by a parametric oscillators. These results hold 
in both cases of single and bichromatic driving of the cavity. In this latter case, 
entanglement becomes larger and more robust against temperature under a 
particular parameter condition in which one mode is driven by a red-detuned 
laser and the other one by a blue-detuned laser. In fact, for equal optome- 
chanical couplings and opposite detunings the system is always stable, even 
for large values of the intracavity power, and entanglement can persist also at 
higher temperatures. 

Finally we have investigated a possible route for coupling optomechanical de- 
vices with atomic ensembles, by showing that if the atoms are placed inside 
the optical cavity and tuned into resonance with the Stokes sideband, optome- 
chanical entanglement is optimally distributed also to the atomic ensemble 
|7U] . Under these conditions one realizes a strongly coupled system showing 
robust tripartite entanglement which can be exploited for the realization of 
CV quantum interfaces [71] . 
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